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ABSTRACT 


A  Steiner  Problem  in  graphs  is  the  problem  of  finding  a  set  of  edges  (arcs)  with 
minimum  total  weight  which  connects  a  given  set  of  nodes  in  an  edge- weighted  graph 
(directed  or  undirected).  This  paper  develops  models  for  the  directed  Steiner  tree 
problem  on  graphs.  New  and  old  models  are  examined  in  terms  of  their  amenability 
to  solution  schemes  based  on  Lagrangean  relaxation.  As  a  result,  three  algorithms 
are  presented  and  their  performance  comparred  on  a  number  of  problems  originally 
tested  by  Beasley  (1984,  1987)  in  the  case  of  undirected  graphs. 


RESUME 

Etant  donne  un  graphe  G  -  (1  ..4)  el  un  sous-ensemble  de  sommets  \ T  de  \  .  le 
probleme  de  Steiner  consiste  a  determiner  un  graphe  partiel  de  G  de  longueur  mini 
male  permettant  de  relier  entre  eux  tous  les  sommets  de  Vl  en  utilisant  event uellement 
un  ou  plusieurs  sommets  de  V  ,1'].  Le  present  article  traite  de  ce  probleme  dans  le  cas 
des  graphes  non-orienles.  Differentes  modelisations  y  sont  examinees,  ainsi  que  des 
methodes  de  solution  basees  sur  la  technique  de  la  relaxation  Lagrangienne.  En  guise 
de  resultats.  tTois  al^orithmes  sont  presentes  et  testes  sur  un  ensemble  de  problemes 
originalement  utilises  par  Beasleyl  1984.  1987  )  dans  le  cas  des  graphes  non-opentes. 
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1.  INTRODUCTION 


1.1  Graph  terminology 

A  directed  graph  G  —  { \ .4}  consists  of  a  tinite  set  of  vertices  1  ( 11  ’;  =  n )  and  a 
set  of  arcs  .4  (  A  =  m  ).  Each  arc  in  .4  is  an  ordered  pair  (r,, t'j)  of  distinct  vertices 
v,.Vj  (.4  C  V  >  \  ).  A  weighted  graph  G  =  { \ ’.  .4 ; ^ .  is  the  graph  G  together  with 
a  nonnegative  real  function  g  defined  over  .4.  Since  our  concern  is  only  with  weighted 
graphs,  we  drop  the  function  g  from  our  notation  for  graphs  unless  otherwise  specified. 

A  subgraph  G1  C  G  is  a  graph  G'  -  {U'-A'}.  where  V'  C  V  and  A'  C  A  and 
such  that  each  arc  a of  G'  {a1^  f  .4* )  is  incident  to  the  same  vertices  in  G1  and  in 
G.  The  weight  of  a  subgraph  G' .  denoted  by  H  (G1 ).  is  the  sum  of  the  weights  cf.  (or 
c,j  for  a j  =  (c(.T  j))  over  all  the  arcs  a in  A1. 

A  *t  mipalh  joining  vertices  rs  and  </  in  G  is  a  sequence  of  vertices  rs.  rs+  iy 
such  that  the  art  (ts_(.  I  f  .4  or  ( i  t_,_  j.  )  €  A  for  each  i.  A  set  of  vertices 

f  1  \  is  weakly  connected  if  there  is  a  semipath  of  vertices  in  V  joining  any  pair 
cf  vertices  in  l  .  A  paih  from  rs  to  iy  in  G  is  a  sequence  of  arcs  in  A  such  that 

(rs.  vs~  j ).  (r<_  ].  t  s_2) . ( J's-i- 1  1 . (vl-l-  V1  )•  The  path  is  simple  if  all  the 

vertices  . r<}  are  distinct.  A  path  is  a  cycle  if  vs  -  iy .  A  set  of  vertices 

5  C  V  is  strongly  connecltd  if  there  is  a  path  from  any  vertex  in  5’  to  any  other  vertex 
in  5  which  does  not  contain  any  vertices  in  \ '  \  S.  A  maximal  strongly  connected  set 
is  a  strongly  conmcied  component  of  G.  An  arhorescence  B,  (B  C  (?)  is  a  set  of  arcs 
such  that 

(i  )  if  ( v, .  tc,  I.  ( ij.  Wj  I  are  distinct  arcs  of  B.  then  u  ,  *  U  j  ;  (al  most  one  arc 
entering  a  vertex  I 

(ii)  B  doe?  not  contain  a  cycle: 

(iii)  B  is  weakly  connected.  {B  is  a  directed  tree.) 

A  graph  G  =  {  \  .  A}  is  said  to  be  connected  if  the  set  \  is  weakly  connected. 
Let  G'  —  be  a  subgraph  of  G  =  {’»  ..4}  and  \  "  C  \  \  then  G'  is  said  to  be  a 

connected  subgraph  with  respect  1c  \  "  if  there  exists  a  semipath  in  G'  between  each 
pair  of  distinct  vertices  in  \  .  Furthermore,  if  G'  is  an  aborescence,  then  subgraph 

G'  is  said  to  be  an  aborescence  subgraph  of  6’  with  respect  to  \  " .  and  if  V1  =  V"  =  1 
then  G'  is  a  spanning  arborescence  of  G.  Note  that  one  verlex  in  V’*  is  designated  as 
the  root  for  the  aborescence. 

1.2  Problem  Statement 

Let  G  =  {1  A}  be  a  weighted  connected  graph.  Let  \  j  C  \'  and  (7(lj  )  =  {6',  C 
G:  such  that  G,  is  a  connected  subgraph  with  respecl  to  l  j}.  The  problem  is  to  find 
the  least  weight  graph  in  G(  \  j ).  Such  graph  (denoted  by  G'*(  \  j  ))  is  called  a  directed 
Steiner  graph  of  G  with  respici  to  l  j .  If  (?*(  1  j  )  is  an  aborescence  then  it  is  called  a 


directed  Sinner  tree  of  G  with  respect  to  V]  or  the  minimal  arborescence  of  G  with 
respect  to  l  ] .  The  vertices  of  the  graph  G’(  l  \  )  r.ot  from  the  set  1  \ .  are  referred  to 
as  Steiner  points. 

Our  graph  notation  attempts  to  tie  together  the  notational  conventions  developed 
by  Bondy  and  Murty  (1976),  Hakimi  (1971)  and  Tarjan  (1977).  An  arc  from  to 
is  denoted  either  by  (r,.t 'j)  or  simply  by  (i,j)  if  no  ambiguity  results. 

We  restrict  the  discussion  to  weighted  graphs  with  strictly  positive  weight  func¬ 
tions  and  to  the  problem  of  finding  the  minimal  arborescence  of  G  with  respect  to 
1 1 ,  i.e..  the  directed  Steiner  tree  problem.  Any  solution  to  this  problem,  designates 
one  vertex  in  \’j  as  the  root  vertex  of  the  resulting  minimal  arborescence  of  G  wilh 
respect  to  lj.  To  facilitate  the  forthcomming  mathematical  formulations,  we  aug 
ment  the  graph  G  with  an  additional  vertex  denoted  as  vertex  0  and  with  a  set  of 

arcs  (O.t,).  for  v,  £  lj.  The  weight  function  g  assigns  to  these  arcs  high  positive 

identical  weights  to  assure  only  one  such  arc  in  the  optimal  solution  to  the  directed 
Steiner  tree  problem. 

The  complexity  status  of  the  Steiner  tree  problem  on  graphs  is  well  settled.  Karp. 
1972.  proves  the  NP-completeness  of  this  problem  by  polynomial  transformation  from 
Exact  Cover  by  3-Sets  problem.  The  problem  remains  NP-complete  if  all  the  arc- 
weights  are  equal,  if  G  is  a  bipartite  graph  with  no  arcs  joining  two  vertices  in  \  \ 

(or  two  vertices  in  \  \  l)  ).  and  also  in  case  G  is  a  planar  graph.  For  more  detailed 

reference  list  on  the  complexity  of  this  problem  see  the  classic  book  by  Garey  and 
Johnson.  1979.  As  far  as  a  literalure  review  of  past  research  on  the  Steiner  tree 
problem  on  graphs,  the  readers  are  directed  1o  the  very  fine  recent  paper  by  Winter. 
1987. 


2.  MATHEMATICAL  FORMULATIONS  WITH  FLOW  VARI¬ 
ABLES 

In  the  mathematical  formulations  that  follow,  the  decision  variables  are  of  two 
basic  types;  (i)  a  binary  (selection)  variables  -  ytJ ,  which  accept  value  1  if  the  aTc 
(i  ,,rj)  is  selected  in  the  solution  and  the  value  0  if  not.  and  (ii)  flow  variables  either 
Tt]r  (or  xtJ)  representing  the  amount  of  flow  through  the  arc  (t>j.r y)  directed  from 
the  root  vertex  0  to  vertex  p  in  1  j  (or  the  total  flow  on  the  arc  (vt.ij)). 

2.1  Mathematical  Formulation  -  A 

Minimize  r  =  ^(l  J )€„  ct]y,}.  0) 

subject  to: 


*  « 


] .  for  f  -  0. 


-jcr_{o}-rijr  -  o}-rj‘r  =  \  0  for '  -  >  for  a11  r  '  1 

- 1 .  for  i  —  p.  -J 


t>jp  $  y>j- 

for  all  (r,.  i -j ) 

€  -4.pt  \\ 

(3) 

T‘JP  z  o. 

for  all  (r,. 

)  €  -4.pt  1  j 

(4) 

y,j  -  0  or  1 . 

for  all  ( i',,  Vj ) 

€  -4 

(5) 

Given  the  objective  of  minimizing  the  total  weight  of  the  selected  arcs  (the  arcs 
for  which  y,j  =  1).  the  constraints  have  to  ensure  the  arborescence  structure  (with 
respect  to  l  j  ).  Constraints  (2)  impose  the  conservation  of  flow  on  all  vertices  in 
V  \('i  -  {0})  and  ensure  that  one  unit  of  flow  leaves  the  root  vertex  0  for  each 
destination  vertex  p  in  \  j.  Constraints  (2)  also  state  that  one  unit  of  flow  reaches 
ea(h  vertex  in  1  j.  Thus,  there  has  to  be  a  path  (in  G)  from  the  vertex  0  to  each  one 
of  the  vertices  in  lj.  Constraints  (3)  ensure  that  flow  is  allowed  only  through  the 
arcs  selected  in  ihe  solution.  Constraints  (4)  are  the  flow  nonnegativity  constraints 
and  (5)  are  the  binary  value  constraints  for  the  appropriate  decision  variables. 

.4  noli :  In  case  the  weight  function  g  of  the  graph  G  —  {V.  .4:  p}  is  not  strictly 
positive,  then  in  order  to  ensure  a  tree  structure  for  ihe  solulion  to  the  MDSTP  we 
add  the  constrains  ytJ  =  1  for  all  j  f  1  j  (2a|.  Since  in  this  paper  we  rest  rid 

the  discussion  to  strictly  positive  functions  g.  constraints  (2a)  are  rtdundenl  and 
subsequently  dropped. 

A  more  compact  mathematical  formulation  to  the  one  presented  in  (A)  is  ob¬ 
tained  by  simply  aggregating  constraints  (3)  with  respect  to  the  index  p.  In  this  case 
we  obtain: 

Lp-:V,  *ijp  S  '  j  y,j-  for  all  €  -4  (3‘) 

We  denote  the  mathematical  formulation  which  contains  the  equations  (1).  (2). 
(3‘),  (4),  and  (5)  as  (A’). 

It  is  of  interest  to  examine  the  implications  of  aggregating  constraints  (3)  into  the 
form  of  (3’).  Since  both  formulations  construct  the  same  weight  solution  if  solved  op¬ 
timal}-,  what  could  be  the  advantage  of  attempting  1o  solve  one  set  of  equations  versus 
the  other  set  of  equations  ?  The  answer  to  this  question  is  in  most  cases  experimental 
in  the  sense  that  a  number  of  researchers  (Magnanti  and  Wong,  1981,  Cornuejols  et 
al.,  1977)  were  much  more  successful!  in  constructing  efficient  solution  schemes  for 
the  disaggregated  formulation  (A)  than  the  more  compact  (A’)  formulation.  In 
order  to  formalize  this  experimental  experience  we  examine  the  linear  programming 
relaxations  (constraints  (5)  become  >  0.  and  ytJ  <  1  for  all  (r,,jj)  G  A )  of  both 
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foriiiuldt ions.  Denote  1  lie  objective  function  value  for  tin-  LI'  relaxation  for  (A)  as 
and  for  (A’)  as 


Theorem  1  :  ; 


_ A ' 
-LP- 


Proof  :  It  is  clear  that  any  feasible  solution  for  the  linear  programming  relax 
ation  of  (A)  is  also  feasible  for  the  linear  programming  relaxation  of  (A’).  In  the 
reverse  direction  1  his  is  not  true  in  general.  One  can  illustrate  it  by  examining  a 
equilateral  triangle  graph  with  4  vertices,  (one  in  the  cenler)  and  7  directed  arcs  (see 
Figure  1  ).  The  l  i  set  consists  of  vertices  1  and  2. 


Figure  1  :  Lqutialerai  triangle  graph. 

The  flow  variables  a"03 1  -  2.  ^3  j  1  =  2.  -r  1 2 1  —  1--T201  -  1  are  clearly  feasible  in 
the  LP  relaxation  for  (A’)  and  not  feasible  for  the  LP  relaxation  for  (A)  • 

Given  two  (different )  equivalent  integer  programming  formulations  (Pi )  and  (P2) 
for  an  optimization  (minimization)  problem.  Denote  by  (RP1)  and  (RP2)  the  respec 
tive  formulations  obtained  by  relaxing  the  integrality  constraints.  If  zpp\  >  ~pP2 
and  not  equal  in  general,  then  formulation  (Pi)  is  called  a  strong  formulation 
with  respect  to  (P2).  This  concept  of  strong  formulation  is  stated  in  Geoffrion,  1979 
and  Van  Roy,  198b. 

In  our  case,  we  proved  that  (A)  is  a  strong  formulation  with  respect  to  (A’)- 

The  next  mathematical  formulation  for  the  minimal  arborescence  problem  with 
respect  to  V’j  requires  only  two  indices  for  its  flow  variables.  It  is  based  on  a  similar 
TSP  formulation  from  Gavish  and  Graves,  1982. 

Mathematical  Formulation  -  B 

Minimize  ;  =  cl}yl} ,  (6) 
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Subject  to: 


-  s: U'.i 

=  - 1 . 

for  all  Ac  l  j 

(T ) 

=  0. 

for  all  A  c  1  \  1, 

(8) 

VI 

for  all  ( 

U.fj)  €  A 

(9) 

iij  >  o. 

for  all  ( 

vi'  vj )  €  A 

(10] 

yt]  =  0  or  1 

for  all  ( 

>,,r;)  £  A 

(11 

The  formulations  (A)  and  (B)  are  equivalent  in  terms  of  determining  the  same 
optimal  solution.  In  terms  of  the  linear  programming  relaxation,  formulation  (A) 
is  a  strong  formulation  with  respect  to  (B)  (the  proof  of  Theorem  1  holds  for  this 
case  too).  What  is  also  clear  is  that  any  feasible  solution  (in  terms  of  flow)  of  the 
LP  relaxation  for  (A’)  is  also  feasible  for  the  LP  relaxation  of  (B)  and  visa  versa. 
Thus.  (A’)  and  (B)  are  fully  equivalent. 

3.  MATHEMATICAL  FORMULATIONS  WITHOUT  FLOW  VARI¬ 
ABLES 

3.1  A  Note  on  the  Minimal  Spanning  Arborescence  Problem  -  (MSAP) 

When  attempting  to  solve  the  MDSTP  on  a  graph,  one  usually  reexamines  the 
very  similar  'easy'  problem  of  constructing  a  minimal  cosl  spanning  arborescence 
(MSAP).  In  all  the  formulations  of  .MSAP  which  are  known  to  the  authors,  flow 
variables  are  used  ( see  for  example  (lavish.  1982)  in  the  same  fashion  as  in  the  MDSTP 
formulations  presented  above  in  Section  2.1.  Aneja,  1980,  presents  a  formulation  of 
the  MDSTP  without  the  flow  variables.  His  is  a  set  covering  formulation  in  which 
the  number  of  constraints  grow>  exponentially  with  the  size  of  problem  instames. 
We  present  a  very  simple  formulation  of  the  MSAP  which  does  not  require  flow 
variables  nor  does  the  number  of  constraints  grow  exponentially  with  the  size  of 
problem  instances.  In  this  new  MSAP  formulation  we  ‘borrow’  a  version  of  TSP 
subtour  elimination  constraints  (Miller.  Tucker,  and  Zemlin,  1960). 

Minimize  r  = 

Subject  to 

Vij  =  1  for  all  >=  1,2,...,  n  (12) 

v,  -  v}  nyl}  <  (n  -  1 )  for  all  i,  j  £  1  ’  U  {0}  (13) 

yl}  -  0  or  1  for  all  i,  j  €  V’  U  {0}  (14) 

where  u,  and  Vj  are  arbitrary  real  numbers. 
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Aote  that  those  cycle  (subtour)  elimination  constraints  (14)  are  for  all  i.j 
l  s  {0}  including  the  'artificial’  root  vertex  0  (contrary  to  i.j  ~  0  in  the  1 SP 
iorinulation )  Also  note  that  the  TSP  subtour  elimination  constraints  of  the  type 
5  ^  y,j  £  S  -  1  and  2  1  for  ab  5  C  V  would  not  be  appro 

priate  for  this  MSAP  formulation. 


3.2  Set  Covering  Formulation  -  C 

First,  we  present  a  set  covering  type  formulation  for  the  MDSTP.  modified  for 
the  directed  graph  by  Wong.  1984.  and  originally  presented  by  Aneja.  1980.  for  the 
undirected  graph. 

Minimize  r  =  A  c‘jV'J  <15) 

Subject  to: 

—  (  v  .  y  y,j  >  1  for  all  .Vj  C  \'  such  that  0  -  A’i  and  .V  j  ^  lj  ^  P 

(16) 


(A  i  denotes  the  complement  of  A’j  in  \  ) 

y,j  fl  or  1  for  all  (i.  j )  -  .4  ( 1 7 1 

In  his  solution  scheme  for  formulation  C.  Aneja.  19^0.  solves  optimaly  the  linear 
programming  relaxation  of  C.  An  interesting  result  due  to  Wong.  19M.  is  that 
:LP  ~  :<LP^  ~IP  - fp b  li  e.,  the  linear  programming  relaxations  of  the  set 
covering  formulation  C  and  the  flow  formulation  A.  have  ihe  same  optimal  values,  i 


The  number  of  constraints  in  formulation  C  is  exponential  in  the  size  of  the 
problem.  There  is  too  little  'structure'  in  this  formulation  to  be  usefull  in  Lagrangean 
relaxation  schemes.  In  order  to  amend  this  'weakeness'  we  add  a  number  of  redundent 
structural  constraints  in  the  next  formulation. 


3.3  Modified  Set  Covering  Formulation  -  D 


Minimize  c  = 

(18) 

Subject  to: 

EIU  y„  -  1 

for  all  j  e  I  j 

(19) 

— (i,j)e  >4  VtJ  -  'c  - 

1 

(20) 

v,  -  Vj  ~  ny,j  <  (n 

-  1 )  for  all  (i.j)  e  A 

(21) 

-  1 

for  all  pel)  and  all  cuts  {('(p)} 

(22) 

where  C(p)  is  a  cut  (a  subset  of  arcs  in  A)  between  the  vertex  0  and  the  vertex 

pel) 


y,j  -  0  or  1  for  all  I i.j)  €  .4 

u,  and  i,j  are  arbitrary  real  number?. 


(23! 


The  objeetive  function  (]S|  and  the  constraints  (22).  (23)  define  the  MDSTP  in 
a  formulation  equivalent  to  C.  Constraints  (19).  (20).  and  (21)  are  redundent.  The 
constraints  in  (10)  ensure  that  only  one  arc  enters  each  vertex  in  \  j .  Constraints  (20) 
state  that  at  least  -  1  arcs  are  selected  in  any  solution  where  initially  \'f  =  V]  and 
in  principle  \]  •_  1 ,  .  Constraints  (21)  are  the  subtour  elimination  constraints.  Con¬ 
straints  (21)  are  not  the  ‘complicating’  constraints  in  this  case.  The  complicating 
constraints  are  the  set  covering  constraints  (22)  which  ensure  a  tree  structure  solution 
which  spans  all  the  vertices  in  \  j .  Dropping  the  constraints  (22)  results,  through  the 
(21  )  constraints  in  a  minimal  weight  forest  solution  to  the  remaining  problem. 

Corollary  1:  :[p  >  :<Lp(  z£p). 

Thi  s  result  simply  follows  from  the  fact  lhal  we  oave  added  constraints  (19)  to 
an  equivalent  formulation  to  C.  The  other  constraints  types  ((20)  and  (21))  would 
be  satisfied  in  the  C  formulation  through  constraints  (16). 

4.  LAGRANGEAN  RELAXATION 

In  this  section  we  describe  a  number  of  Lagrangean  relaxations  for  the  mathe¬ 
matical  formulations  presented  in  Sections  2  and  3.  Lagrangean  relaxation  approach 
for  solving  ’hard’  problems  js  based  on  the  observation  that  by  removing  the  compli¬ 
cating  constraints  from  a  mathematical  formulation,  the  resulting  problem  is  ‘easily’ 
solvable.  A  solution  to  the  relaxed  problem  constitutes  a  lower  bound  on  the  solution 
to  the  original  problem.  The  thrust  in  such  an  approach  is  to  obtain  a  maximal  lower 
bound  which,  if  it  does  not  solve  the  original  problem,  can  be  integrated  into  an 
implicit  enumeration  scheme  such  as  branch  and  bound. 

4.1  Lagrangean  Relaxations  of  (A) 

We  present  two  Lagrangean  relaxations  of  formulation  A.  In  the  first  one  the  arc 
selection  constraints  are  relaxed,  resulting  in  a  shortest  path  type  problem.  In  the 
second  relaxation  the  conservation  of  flow  constraints  are  moved  into  the  objective 
leading  to  a  more  difficult  subproblem. 

4.1.1  The  First  Relaxation 

In  formulation  (A),  the  complicating  constraints  are  the  arc  selection  constraints 
(3).  which  ensure  that  the  unit  flow  from  the  root  vertex  0  to  a  vertex  p, p  <=  Ij  passes 
only  through  the  arcs  selected  in  the  solution.  Below  we  present  the  Lagrangean 
relaxation  obtained  by  moving  the  constraints  (3)  multiplied  by  nonnegative  A XJp 
into  the  objective  function.  For  a  given  veclor  of  multipliers  A.  the  problem  is 

-  Q  ■ 


Minimize  -(A)  —  —  (i,j)€-4  — 1  i.ji)€.4  ^-pC)  i  ^ijp^Vij  X>]P^ 

which  after  rearrangement  of  terms  has  the  following  form 

Minimize  . ( A )  —  X]( i.j)£ .4 ( cij  ~  1 1  'Vp  )ibj  ~  A  ^/>£ 1 1  jpx  >jp  ( 24 ) 

Subject  to: 


el.  for  l  =  0,  *N 

T  -  V'  .. 

■VP  — Jtl 

o{0}  XJIP  -  \  0-  for  J  4  0,p  >  for  all  p  £  \  j 

(25) 

v  1 .  for  i  —  p,  J 

VI 

ft, 

H 

V'l 

o 

for  all  (i,j)  €  -4,p  £  Ij 

(26) 

y,j  =  0  or  1. 

for  all  (i.j)  £  A 

(27) 

The  'best'  Lagrangean  value  is  obtained  by  maximizing  c(A)  over  nonnegative 
A 's  (A ,jj,  >  0.  for  all  (i.j)  £  A.p  £  l)).  For  those  optimal  Lagrangean  multipliers 
(A jjj,)  the  following  relationship  holds  (see  Gavish,  1978.  and  Wong.  1984.  for  the 
dual  formulation  of  A): 

G?  “  UpeG  Vr  >  0  for  all  (i.j)  £  .4  (28) 

By  observing  that  the  above  ((24)  ■  (27))  Lagrangean  formulalion  has  ihe  Inte¬ 
grality  Property  (Geoffrion.  1974).  the  lower  bound  value  oblained  for  the  MDSTP  by 
solving  (24)  -  (27)  is  equal  to  the  value  for  the  linear  programming  relaxation  to  the 
problem.  In  this  case,  the  main  advantage  for  examining  the  Lagrangean  relaxation 
of  A  would  depend  on  how  fast,  in  comparison,  can  such  a  relaxation  be  solved.  In 
addition,  such  a  solution  could  be  more  amenable  for  developing  good  heuristics. 

The  ( 28 )  inequalities  suggest  a  fast  solution  procedure  to  obtain  the  maximal  ;(  A  j 
solution  for  the  Lagrangean  relaxation  of  A  based  on  the  repeated  use  of  the  shortest 
path  algorithm  with  cost  modifications  along  the  way.  In  the  algorithm  outlined 
below,  we  succssively  adjust  the  values  of  the  A, ^"s  using  the  subgradient  method 
described  in  Held  et  al.  (1974)  while  preserving  the  dual  feasibility  of  these  multipliers 
via  the  (28)  inequalities.  We  enforce  for  all  (i,j)  £  A  throughout  the  algorithm 
below  (in  each  modification  of  multipliers)  the  following  constraints  :  Vp  ~ 

ct].  Then  the  objective  of  the  Lagrangean  problem  is  equivalent  to  a  shortest  path 
problem. 

Denote  by  z(LAl(p))  the  value  of  :{LA\)  obtained  by  the  Algorithm  LAI  outlined 
below,  given  that  p,  ( p  £  )  is  the  root  node  of  the  Steiner  tree.  Denote  by  :(LAl(p)) 

the  weight  of  the  tree  generated  by  :{LAl(p))  solution  (i.e..  assign  the  actual  costs 
c,j  to  the  arcs  in  the  Steiner  tree).  Let  Sp{l.p )  denote  the  shorlest  path  from  /  to  j< 
in  the  network  with  arc  costs  (A,^) .(i.j)  £  A. 


■  10  - 


Algorithm  LAI 

Step  0  :  (Initialization) 

Set  \l}},  -  for  all  (i.j\  *£  A.p  €  V\ 

:(^)best  =  O'^bist  =  oc:  I,er  -  0.  ImpJier  -  0; 

Idivide  -  (preset  parameter);  b  =  2;  Iterlim  =  Maximum  number  of  iterations; 
tf,  —  smaller  gap  on  <?;  t£,es<  =  relative  precision  desired  on  ;  st  =  0 

Step  1  :  Solve  the  following  (Lagrangean  relaxation)  problem: 

min  Z]  =  Ylp£\ H(,  j)£A  rx*ir  sut>ject  to  constraints  (26)  and  (27). 

This  problem  can  be  solved  with  a  shortest  path  algorithm  as  follows: 

For  each  p  6  l  j  compute  all  the  shortest  paths  from  p  to  k  E  Vj  using  arc  costs 
{A,,,,}  and  denote  its  cost  by  S},.  Selecl  the  solution  with  the  minimal  Sp  value. 

Step  2  :  U  pdating  the  bounds 

Lower  bound  : 


If  z{\)  >  :l  A  o  then 

-(A  )(l(sl  ^  :(A).  and  for  all(i.j)  i  A.p£  Vj.A  ,Jp(best)  -  \,jr.~i,jr(besi )  = 
st(be  si )  =  st .  Implter  -  0 

Upper  bound  : 

Consider  the  subgraph  C  =  (U  _  { 0 } .  ,4 r )  where  A'  =  {(i . ]  )■  x,JP  >  0  for  some 
p  £  Uj }.  Obtain  the  shortest  path  tree  structure  for  this  subgraph  by  computing  the 
shortest  paths  from  0  to  each  node  of  l  j.  Let  { =  1}  for  each  arc  ( i.j )  in  this 
shortest  path  tree  and  0  otherwise. 

The  current  solution  vector  (y,x)  is  always  feasible  for  the  original  problem  and 
its  cost  Zf{y,r)  =  OjVij 

If  -f(y-T)  <  -besl  lhen 

-best  =  zf(y-*),  Implter  =  0 
Step  3:  Updating  A 

(i)  If  Implter  =  Idivide  then 

Change  the  value  b  and  restart  from  the  best  solution  c(A )(K$t  (i.e..  b  =  b/ 2, 
Implter  =  0?  r(A)  =  z(X)hes1.st  =  st(best)i2.  \,}p  =  Aup(6esf  JOijp  =  hjp(b^st)  for 
all  (i,j)  €  A.p  €  Uj ). 


■  u  - 


Eht  Compute  the  ascatt  dtrtchou  y  and  the  step  si  for  the  current 

solution. 

T  *,)i  ~  y,j-  for  all(i.j)  £  A.p  ‘z  \  i 
s I  _  ^(-S;»f~-(A)) 

A  " 

( it )  ^ijji  —  ~  !)]••  ( 1  •  ]  I  £  -4.pt  1  ] 

The  new  multipliers  \,jj,  are  obtained  by  solving  the  following  problem: 
min{  A  -  A  =  ciy\lJT,  >  0.(t,j)  €  -4.p  £  Vj } 

The  vector  A  is  the  projection  of  A  on  A,);i  =  c,rAlj;,  >  O.(i.j)  e  A.p  £ 

ll) 

.\of(  ■  The  solution  of  the  projection  problem  follows  the  procedure  suggested 
by  Held  et  al.  (197-1.  pp.  771. 

Si <- p  J,  :  Sloping  Conditions 

(a  I  if  lit  r  -  JUrlim  then  Stop 

( b  l  If  <.  t  hen  Stop 

(rl  lf  <  ,hen  S,G{’ 

Otherwise  Go  To  Step  2 

4.1.2  The  Second  Relaxation 

Following  a  different  relaxation  approach,  we  dualize  on  the  flow  constraints  (2) 
using  A ,j,  as  the  Lagrangean  multipliers.  The  number  of  multipliers  reduces  to  7?  I  j  . 
which  is  considerably  less  than  in  the  previous  relaxation.  The  objective  function 
becomes: 

Minimize  r(A)  =  £(l,»e.4  cijV»j  -  *  0}P  -  T.jO'rj0p  ~  1 )  - 

— P<E '  i  ^Pp(  1  Tpjp~^L  j£  f  r jpp +  ^  )—  —  ig  V.«^0,p  ^-p£  1 1  ^»p(  V  T)tp ) 

After  rearrangement  of  terms  the  objective  is  converted  to: 

Minimize  -(A)  —  ci jVij  ~+‘  -3p€' \  — (i j)€ .4  •GjpGjp  ~  ^-p£l’i  ^0 PP  (30) 

Subject  1o  (3),  (4),  and  (5)  where  cIJ},  -  \}},  -  for  all  (i.j)  £  .4  and  p  £  1 ) . 

Note  that  the  A^'s  in  this  relaxation  are  unrestricted  in  sign  and  the  can  be 
handled  implicitly  storing  n‘  entries  instead  of  r?3  entries. 


In  order  to  strenghten  the  lower  bound  obtained  from  such  a  relaxation  we  amend 
it  with  the  following  constraints: 

—  D>}  -  1  for  all  J  v  '  l  (30) 

and  x,jj  =  y,j  for  all  \i.j)  -  A  and  j  t  l)  (31) 

The  constraints  (30).  (31)  are  redundent  in  the  formulation  A,  but  are  helpful 
in  increasing  the  Lagrangean  bound.  Note  also  that  (30)  is  a  relaxation  of  ihe  (tree) 
constraint  which  ensures  that  only  one  arc  enters  a  node  in  V]  (i.e.,  ^i£l'  Vij  —  1  f°r 
all  j  c  Vj ).  The  'tree'  constraint  is  tighter  but  nol  easily  solvable. 

Before  presenting  a  solution  procedure  for  the  problem  defined  by  (29).  (3),  (4). 
(5).  (30).  and  (31),  we  make  the  following  observations: 

(i )  From  the  selection  constraints  (3)  we  notice  that: 

(a)  If  y,j  -  0  then  =  0  for  all  />  £  1 1 

(b)  If  y,j  =  1  then 

(I)  i,jj  -  1  which  follows  from  (31).  and 

(II  I  tIJ}i  =  1  if  rl]},  •-  0.  and  =  0.  if  c., J;,  2  0  A  /»  ~  J- 

Note  that  for  the  same  reason  as  in  the  first  relaxation  of  A.  the  maximum  lower 
bound  for  the  MDSTP  oblained  by  solving  this  relaxation  can  nol  exceed  ihe  bound 
obtained  from  the  linear  programming  relaxation  of  the  problem. 

Algorithm  LA2 

Sit}'  1:  (Initialization  I  Set  A ,p  =  :  Iter  =  0:  Implter  =  0:  ;( A =  0:r{,(S/  = 

cosl  (>f  the  best  feasible  solution  found  so  far:  Idivide  =  presel  parameter:  t  -  2:  si  - 
0:  Ilerlim  -  maximum  number  of  iterations:  tf,  =  smaller  gap  on  -  relative 

precision  desired  on  -( ^  )b(st- 

Note  that  A®  are  randomly  generated.  is  oblained  by  computing  the  shortest 
paths  from  one  node  of  1 )  to  all  other  nodes  in  1 j . 

Step  2:  Solve  the  Lagrangean  Problem: 

Iter  =  Iter  -f  1.  Implter  =  Implter  -t-  1.  Set  yt]  —  0.x, jp  =  0.  for  all  (i.j)  € 
-4 ,p€  V). 

(1)  For  each  arc  (i.j)  €  .4 

(i)  Compute  M,j  =  ctJ  +  min{0:cIj;)} 

(ii)  If  j  €  V]  set  M,j  =  MtJ  -  cl}] 

(iii)  If  M,j  <  0  then  set  ytJ  =  1  and 


•  IS  - 


For  all  p  ■_  Vj  set 


^.1  for  clJp  <  0  or  j  -  p 

T>JV  =  < 

v0  otherwise 

(2)  For  each  j  £  Vj.  if  MtJ  >  0  for  all  (i,j)  €  .4  then  determine  /  such  that 
Mjj  =  min,.f( }  and  set  yij  =  1.  for  all  p  such  that  qJp  <  0  or  p  =  j  set 
27,.,  -  1.  otherwise  set  xt]p  =  0. 

(3)  Compute  -(A) 

Si (p  S:  (l  pdating  the  bounds) 

( 1  )  If  ct  A  1  ;■  ;(  then 

-I  ^  =  r(A).At(S<  =  A.  ~i,p(l'(s1 )  =  -;,p.  St  (best )  -  st.  Impiter  =  0. 

(2)  If  the  solution  (y.r)  is  feasible  for  the  original  problem  then  compule  its  cosl 
and  denote  the  value  by  zp. 

If  Zf  •  Z(l(t1  then  set  zi,tSt  ~  zf.  Set  Impiter  =  0: 

St<p  4:  (l  pdating  the  A) 

(i)  If  Impiter  -  Idivide  Ihm  set  A  =  Impller  =  0.  z(  A)  =  zl\)i„s1.st  - 

sl(htst)  2.  A  ip  -  si ).  ~),j  -  ->(pl  be  sf 1  for  all  i  ~  l./>€  \] . 

fist  compule  the  ascenl  direction  *  and  the  si ej»  st 

A'Lj-lY  *0jr~  — >€l  Tj0p  ~  1  >  f°r  all  p  €  \] 

“tip  i-j^v  *IJP  -  Y.jiV  Tjip)  for  all  /-  -  r,.(  zt  O.p 

(Hjtl  rTJP  ~  — l  ■rjpr>  ^  ^  f°r  P  ^  f  1 

<7  =  /; 

*>; 

(ii )  ^ ip  ~  $1  ~)tp- 1  £  f  •  p  €  f  1 


Step  5:  (Sloping  condition) 

(a I  If  Iter  >  lierlirmi 

Stop 

c r 

-v* 

A 

Stop 

(c) , f 

Stop 
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Otherwise  Go  To  Step  2. 


4.2  Lagrangean  Relaxation  of  (D) 

The  complicating  constraints  in  the  D  formulation  are  the  set  covering  -  (22) 
constraints.  The  major  difficulty  is  that  the  number  of  constraints  in  (22)  is  expo¬ 
nential  in  the  size  of  \  j  set.  On  the  other  hand  the  number  of  ytJ  variables  is  exactly 
n(n  -  1 ) ,  2 .  This  implies  that  in  the  linear  programming  relaxation  of  this  formula¬ 
tion.  most  of  the  constraints  in  (22)  are  nonbinding.  The  difficulty  lies  in  finding  the 
binding  constraints.  By  removing  the  constraints  (22)  and  adding  them  to  the  objec¬ 
tive  function  multiplied  by  the  appropriate  (nonnegative)  Lagrangean  multipliers  we 
obtain  the  following  relaxed  formulation: 

Minimize  -(A)  =  c,}ytJ  ~  HC|l€n  V1  ~  C(/>)  JO; ) 

Subject  to  constraints  (19).  (20).  (21)  and  (23)  where  II  is  the  index  set  of  all 
cuts  Cj,  between  the  vertex  0  and  the  vertices  p  £  lj. 

This  formulation  can  be  rewritten  as: 

Minimize  :(A)  ^  A  ctjy>j  (32) 

Subject  to  constraints  (19i.  (20).  (21).  and  (23)  where: 

»o,  j  -  -V,  for  ii.j)  ('(;>)  for  some  p  -  \  j 


l  c,j  for  i  i.j)  f  ("'(pi 

This  last  problem  ((32).  (19).  (20).  (21).  and  (23))  for  fixed  values  of  A  can  be 
solved  in  polynomial  time.  (Again  from  Lagrangean  duality  we  gel  c,j  >  0  for  all 
(i.  j)  £  .4.)  For  a  given  A  vector  denote  this  problem  by  LDl. 

The  question  now  reduces  to  one  of  finding  A*  which  maximizes  the  value  of  r ( A ) . 
i.e..  r(  A* )  =  A/arA>0{:(A)}. 

Theorem  2  :  r(  A* )  > 

Proof  :  This  result  is  based  on  the  observation  that  the  set  of  constraints  (19). 
(20),  (21).  and  (23)  do  not  have  the  Integrality  Property  (GeofFrion,  1974,). 

4.2.1  Dual  Ascent  Procedure  for  Initial  Ap  Values 

In  this  section  we  describe  a  dual  ascent  procedure  for  computing  a  lower  bound  to 
the  undirected  version  of  the  MDSTP  using  the  D  formulalion  of  the  problem.  Mod- 
ifing  this  dual  ascent  procedure  for  a  general  directed  graph  is  left  as  an  algorithmic 


exercise. 


First  we  explain  and  provide  an  outline  of  the  algorithm  followed  by  a  small 
numerical  example  and  conclude  with  a  detailed  description  of  a  dual  ascent  procedure 
for  finding  good  multipliers  to  LDl. 

The  algorithm  begins  by  solving  the  problem  LDl  with  -  {0}  _  \'j  and  all  A 
values  equal  to  zero.  Let  Lg  be  the  cost  of  this  solution.  The  network  structure  of  this 
solution  is  that  of  a  ‘sparse'  forest  F  of  disconnected  components.  In  ^ase  of  a  directed 
graph,  this  problem  can  be  solved  with  a  modified  Tarjan's  algorithm  (Tarjan.  1977). 
We  differentiate  between  two  component  types.  Components  which  conlain  at  least 

one  node  belonging  to  Vj  {0}  are  denoted  by  T  =  {7q.  T\.  To _ }.  where  Jq  is  the 

component  which  contains  the  root  node.  The  second  set  of  components  consists  of 
the  points  S.  5  C  V  -  \’j  not  contained  in  any  of  the  components  in  T. 

At  this  point  we  pick  a  component  Tk..k  ?  0.  and  compute  the  minimal  cost 
of  expanding  this  component.  Let  />/.  =  Min  cjk  rfy  jk  {c,j } .  6k  is  ihe  multiplier 
value  over  the  cut  seperating  the  nodes  in  T k  from  all  the  other  nodes,  ^he  cost 
matrix  is  updated  1o: 


.c,j  -  for  all  j  e  J\  .i  7  l'  Tk 

c,j  -  < 

*  C)j  otherwise 

As  a  result  of  cost  matrix  update,  one  or  more  arcs  on  ihe  cul  have  a  reduced 
cost  of  zero. 

Let  L  gt  -  Lgt  ,  -  Sk  where  t  is  the  iteration  number. 

The  component  is  merged  via  the  zero  cost  arcs  with  a  number  of  other 
components.  This  process  continues  till  Tk  is  merged  with  Tq  al  which  point  anolher 
component  in  T  is  selected.  At  each  iteration  al  least  one  component  is  added  to  Tk 
thus  this  merging  of  components  stops  at  mosl  after  5  -  T  -  1  steps.  When  the 
process  is  completed  we  have  only  one  component  Tq  which  contains  a  subset  5j  of 
nodes  (Steiner  points).  $’j  C  \  \  Vj .  We  remove  from  Tq  all  the  nodes  6  5j  which 
are  dangling  nodes  (i.e..  their  degree  is  <  1).  Let  5’q  be  the  set  of  terminal  nodes. 

LB  =  LB,  -  Hf£S(,cp,.s 

where  ps  is  the  adjacent  node  of  node  5..s  €  Tq.  This  last  ‘trimming’  step  is 
repeated  until  all  the  terminal  nodes  (degree  <  1)  are  nodes  in  Vj  only. 

lb  is  the  low-er  bound  value  to  the  MDSTP. 

In  order  to  illustrate  those  steps  we  use  the  following  example.  The  example 
consists  of  3  required  nodes  (1.  3.  and  4.  where  node  1  is  designated  as  the  root  node  I 
and  4  potential  Steiner  points.  The  network  matrix  (the  arc  weights)  is  symmetric, 
which  in  the  directed  graph  version  implies  two  arcs  of  the  same  weight  and  oposite 
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direction  between  two  adjacent  nodes  of  the  graph.  The  initial  lower  bound  obtained 
for  this  example  is  5.  The  lower  bound  after  two  multiplier  adjustments  is  8.  while 
the  optimal  solution  to  ’he  problem  is  9. 


Example  1  : 
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Figure  2  :  The  original  'distance' matrix. 

The  initial  'forest'  solution  connects  node  1  to  node  2.  node  3  to  node  5.  and 
node  4  to  node  6.  This  solution  has  a  value  of  5.  We  pick  the  tree  containing  the 
nodes  3  and  b  together  wiih  the  arc  from  3  to  b  as  our  7*  Iree  in  the  algorithm.  The 
corresponding  <l'j  value  is  3  (the  arc  weight  froiti  node  2  to  node  3)  and  the  new  lower 
bound  is  b  -3  =  8.  The  modified  distance  matrix  is  as  follows: 
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Figure  3  :  The  modified  ‘distance’  matrix  after  one  descent. 

A  new  tree  is  selected  since  the  previous  one  contains  the  root  node  1.  The 
tree  is  the  arc  from  node  4  to  node  6  together  with  the  two  nodes.  The  new  fy, 
corresponds  to  the  arc  from  node  3  to  node  4  and  has  a  weight  of  2.  The  new  lower 
bound  is  5  -  3  -  2  =  10.  Since  the  expanded  new  tree  contains  the  root  node  and 
there  are  no  trees  remaining,  we  trim  the  Iree  from  the  dangling  not  required  nodes 
and  obtain  a  tree  which  contains  the  nodes  1.  2.  3,  and  4  and  the  arcs  (1.2),  (2,3). 


and  (3.4)  for  the  total  of  10  -  2  =  8.  This  solution  corresponds  to  a  feasible  Sieiner 
tree  of  value  of  11  instead  of  the  optimal  value  9. 
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Figure  4  :  The  final  modified  'distance'  matrix. 

4.2.2  The  Dual  Ascent  Algorithm 

Step  1  :  We  start  with  all  A, ;  =  0  and  solve  the  problem  ((32).  (19).  (20). 
(21).  and  (23))  to  obtain  z(0).  The  solution  to  this  problem  is  in  the  form  of  a  not 
necessarily  connected  set  of  trees  T  which  might  not  include  all  the  nodes  in  Y  (i.e.. 
a  'sparse'  forest  ).  In  case  we  obtain  only  one  tree  then  we  have  the  oplimal  solulion 
for  the  Sieiner  tree  problem.  Denote  by  L g  the  cost  of  this  solution. 

Slip  J  :  Let  Y  be  the  number  of  trees  in  the  foresl  and  let  7)  be  the  set  of 
nodes  in  tree  k.  One  of  these  trees  conlains  the  root  node.  Denote  that  tree  by  7p. 
Denote  by  S  the  set  of  nodes  in  1  .  )=(|7j  (i.e..  the  potential  new  Steiner  points). 

Step  3  :  Pick  one  of  the  trees  in  T  7U .  Tree  Tm  for  example.  Oompule 
-  =  Minj£Tn.t€Tm{z'ij} 

Let:  Lg  -  Lg  t  L,  c,j  =  c,j  -  L.  for  all  i  €  Tm,J  g  Tm.  or  j  €  Tm.i  g  T„, 

Step  J,  :  (  Merging)  Every  component  in  T  \  Tm  and  S  with  exactly  one  zero 
cost  arc  to  Tm  (i.e..  ctJ  =  0  j  €  Tm,i  g1  Tm )  is  merged  with  Tm  creating  a  new 
forest.  (In  case  of  components  with  multiple  zero  cost  arcs  to  Tm  see  Remark  1.) 
Rename  the  trees  in  the  new  forest.  Repeat  Steps  3  and  4  until  Tm  is  merged  with 

To- 


Step  5  :  If  T|  g  1  then  go  to  Step  2,  otherwise:  eliminate  all  the  terminal 
nodes  in  7q  which  belong  to  V  \  Vj  and  reduce  the  corresponding  Lg  value  by  the 
corresponding  arc  costs.  I.e..  Lg  =  Lg  -  cp,  s-  Tg  is  the  lower  bound  value 

for  the  Steiner  tree  problem. 

Remark  1  :  If  we  add  multiple  zero  cost  arcs  between  Tm  and  another  tree  in 
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T  then  a  cycle  is  crealed.  Thus,  at  each  merge  operation  only  one  arc  can  be  added 
between  any  pair  of  trees.  Since  without  loss  of  generality  we  can  assume  (different  1 
integer  arc  weights,  we  can  expect  thal  the  number  of  zero  cosl  arcs  between  a  pair 
of  trees  in  a  merge  slep  is  small.  This  suggesls  a  parallel  processing  algorithm  for 
the  construction  of  the  new  trees  (merged  trees)  at  Step  5.  In  each  case  were  more 
than  one  (say  k)  zero  cost  arc  exists  between  a  pair  of  trees,  k  new  merged  trees  will 
be  stored  and  processed  in  parallel.  At  the  end.  following  Step  5.  we  obtain  the  best 
lower  bound  by  examining  all  the  trees  grown  in  parallel. 

5.  COMPUTATIONAL  RESULTS 

The  solution  methods  which  were  developed  and  described  in  the  previous  sections 
provide  a  lower  and  upper  bounds  on  the  value  of  the  optima]  Steiner  tree  solution 
for  a  given  graph.  Optimal  solutions  were  obtained  for  the  problems  for  which  the 
difference  between  the  value  of  the  upper  bound  and  the  lower  bound  was  less  than 
one.  In  order  to  investigate  the  comperative  performence  of  the  solution  methods 
developed  in  this  paper,  they  were  programmed  and  tested  on  a  set  of  problems 
taken  from  Beasley  (]9M|.  We  present  the  results  of  these  tests  in  TABLE  I 
below.  The  algorithms  LAI  and  LA2  were  programmed  in  FORTRAN  and  the  Dual 
Ascent  Algorithm  (D.A.  Algorithm  in  Table  I)  was  programmed  in  PASCAL.  The 
data  set  is  the  one  tested  in  Beasley  (19*4)  and  consists  of  18  randomly  generated 
problems  for  undirected  graphs.  For  algorithms  LAI  and  LA2  we  have  considered  the 
directed  version  of  these  problems  by  duplicating  each  arc  and  assigning  directions. 
The  three  algorithms  were  tested  using  a  VAX  8600.  In  all  the  tests  the  number  of 
subgradiant  iterations  was  restricted  to  800.  the  initial  £  value  was  set  to  2  and  the 
parameter  Idivide  was  sel  to  20  for  LAl  and  to  40  for  LA2. 
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Problem 

Number 

D 

A 

V, 

Algorithm  LAI 
Lower  Cpper 

Bound  Bound 

Algorithm  LA2 
Lower  Upper 

Bound  Bound 

D.A.  Algorithm 
Lower  Upper 

Bound  Bound 

1 

bO 

126 

9 

O 

81.44 

82* 

72 

83 

2 

50 

126 

13 

83.00 

83* 

82.88  ; 

90 

60 

124 

3 

50 

126 

25 

137.88 

138* 

137.07  ' 

177 

107 

160 

4 

50 

m 

9 

59.00 

59* 

58.95  | 

77 

48 

93 

5 

50 

BB 

13 

60.98 

61* 

60.58  ! 

65 

49 

109 

6 

50 

iii 

25 

121.65 

122* 

121.33 

148 

84 

154 

/ 

H 

188 

13 

110.97 

111* 

108.79  1 

123 

92 

138 

8 

B 

1 SJ' 

19 

103.99 

104* 

101.81 

118 

73 

125 

9 

188 

38 

219.79 

220* 

207.81 

234 

207 

255 

75 

300 

13 

85.90 

86* 

84.81 

124 

61 

175 

11 

1  r> 

300 

19 

88.00 

88* 

87.79 

127 

64 

192 

12 

To 

300 

38 

172.20 

3  74 

166.20 

228 

139 

204 

13 

250 

17 

165.00 

165* 

162.77 

192 

94 

263 

14 

100 

"ol 

_ 

2  b 

234.91 

235* 

224.56 

277 

131 

310 

lb 

Bt§ 

250 

317.60 

318* 

301.27 

353 

249 

387 

16 

m 

17 

127.00 

127* 

122.65 

162 

73 

256 

17 

MM 

2b 

128.17 

131 

124.52 

143 

101 

163 

18 

im 

400 

50 

215.56 

218 

209.52 

280 

■  - 

182 

260 

TABLE  I  :  Computational  results  for  the  three  algorithms  for  Sleiner  iree 
problem  on  graphs. 

Out  of  the  18  problems  attempted,  lb  were  solved  optimally  by  the  LAI  algorithm. 
Only  the  first  problem  was  solved  optimally  by  the  LA2  algorithm  and  the  Dual 
Ascent  algorithm  did  not  produce  a  single  optimal  solution.  In  terms  of  the  quality 
of  the  lower  bound  values,  LAl's  lower  bound  values  dominate  the  values  generated 
by  LA2  and  the  dual  ascent  algorithm.  (The  optimal  solution  is  noted  by  * . ) 


SAMMERY  AND  CONCLUSIONS 

We  have  presented  a  number  of  mathematical  formulations  for  the  directed  and 
undirected  Steiner  tree  problem  on  graphs.  These  formulations  have  been  used  to 
develop  Lagrangean  based  lower  bounding  procedures  for  the  problem.  In  compu¬ 
tational  tests  (on  18  problems  used  by  Beasley  (1984.  1987)  for  testing  undirected 
Sleiner  tree  problems),  it  has  been  shown  that  one  of  the  algorithms  (LAI )  generates 
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lower  bound  values,  that  are  close  lo  the  optimal  solutions.  The  norifeasible  solutions 
generated  by  the  Lagrangean  based  procedure  have  been  incorporated  into  heuristics 
which  attempt  to  generale  "good"  feasible  solutions.  Here  again  algorithm  LAI  has 
generated  optimal  solutions  to  la  out  of  the  18  problems.  For  the  other  3  problems 
the  gap  between  the  feasible  and  lower  bound  values  were  under  2/i .  Note  that  in 
case  of  an  undirected  graph.  Beasley  ( 1984  )  reports  solving  to  optimality  only  G  out  of 
the  18  problems,  and  Beasley  (1987)  again  for  the  undirected  graphs  reports  solving 
to  optimality  15  out  of  the  1*  problems  using  Cray  X-Mp  48  machine.  The  combined 
results  (Beasley  and  ours)  solve  optimally  17  out  of  the  18  problems. 

Based  on  the  above  results,  it  is  our  believe  that  algorithm  LAI  can  be  used  as 
an  effective  tool  in  Branch  and  Bound  based  procedures  for  solving  the  problem. 


•«  l>  'lynx  nl  :  We  thank  J.E.  Beasley  for  providing  us  with  a  copy  of  his  tesi 

problems. 
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